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The Twin-Jet System in NGC 1052: 
VLBI-Scrutiny of the Obscuring Torus 
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Abstract. NGC 1052 offers the possibility to study the obscuring torus around a supermassive black hole, predicted 
by the standard model of active galactic nuclei, over a wide range of wavelengths from the radio to the X-ray regime. 
We present a detailed VLBI study of the parsec-scale structure of the "twin-jet" system in NGC 1052 in both total 
and polarized intensity and at multiple frequencies. We report the detection of linearly polarized emission from the 
base of the eastern jet at 5 GHz. While the radio spectrum in this region might be still consistent with synchrotron 
self absorption, the highly inverted spectrum of the western jet base represents a clear sign of pronounced free- 
free absorption in a circumnuclear torus. We observe an abrupt change of the brightness temperature gradient 
at a distance of ~ 0.2 pc to 0.3 pc from the central engine. This might provide an observational signature of the 
edge of the central torus, where the transition from an external pressure-dominated jet regime to a more or less 
freely expanding jet takes place. We determine the absorbing column density towards the western jet core to be 
~ 2.2 X 10^^ cm~^ in good agreement with the values derived from various X-ray observations. This suggests that 
the nuclear X-ray emission and the jet emission imaged by VLBI originate on the same scales. 
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1. Introduction 

Radio observations of the low-luminosity active galac- 
tic nucleus (AGN) NGC 1052 with Very Long Baseline 
Interferometry (VLBI) at multiple frequencies have re- 
vealed the presence of a dense circumnuclear absorber, 
which obscures the very center of this elliptical galaxy 
(Kellermann et al. 119991 Kameno et al. I2(JU1I Vermeulen 
et al. l2(JU3fl . Indeed, the standard model of AGNs predicts 
the existence of an obscuring torus, whose inner surface 
is expected to be photo-ionized by illumination from the 
accretion disk. Colder, neutral material forms the outer 
boundary of the obscuring torus. NGC 1052 provides the 
possibility to study the physical properties of this obscur- 
ing torus complementary in various wavelength regimes 
and with a variety of observational methods. Particularly, 
the combination of VLBI and X-ray spectroscopic stud- 
ies is capable of addressing the same basic questions with 
complementary methods. Various X-ray observations of 
NGC 1052 imply a model-dependent column density of 



ray core (Guainazzi et al. 119991 Weaver et al. 119991 Kadler 
et al. 2004 ) but the angular resolution of X-ray telescopes 
is not sufficient to measure the accurate position and ex- 
tent of the absorber. Relativistically broad iron line emis- 
sion at 6.4 keV is seen in a high-quality X-ray spectrum 
obtained with the XMM-Newton telescope (Kadler et al. 
in prep.). Thus, being the first radio-loud AGN with a 
strong compact radio core that exhibits strong relativis- 
tically broadened iron line emission from the inner ac- 
cretion disk, NGC 1052 provides a unique possibihty to 
study the inter-relation between AGN mass accretion and 
jet-formation. For future combined VLBI structural and 
X-ray spectroscopic monitoring observations, it is essen- 
tial to study in detail the influence of the obscuring torus 
on the parsec-scale jet structure at radio wavelengths. 

NGC 1052 is a moderately strong, variable source in 
the radio regime, with a luminosity (integrated between 
1 GHz and 100 GHz) of - 4.4 • lO^'^ergs"! fWrobel lTOMjl 
and an unusually bright and compact radio core. Together 
with the proximity of the source of only 22.6 Mpc^ this 
makes NGC 1052 a premier object for VLBI studies, aim- 
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^ We assume a Hubble constant of Ho — 65kms~^Mpc~^ 
and use the measured redshift of z = 0.0049 by Knapp et al. 
(I1978II . This results to a linear scale of 0.11 pcmas"^. 
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ing at the ultimate goal of revealing the physical properties 
of obscuring tori in AGNs. The parsec-scale structure of 
NGC 1052 shows a twin jet with an emission gap between 
the brighter (approaching) eastern jet and the western (re- 
ceding) jet and free-free absorption towards the western 
jet (Kellermann et al. 11999(1 . Kameno et al. (|2U(J15 sug- 
gested the presence of a geometrically thick plasma torus 
and a geometry of the jet-torus system in which 0.1 pc 
of the eastern jet and 0.7 pc of the western jet are ob- 
scured. The multi-frequency VLBI structure of NGC 1052 
has been further studied by Kameno et al. (2003) and 
Vermeulen et al. H2003fl . The kinematics of both jets at 
2 cm have been investigated by Vermeulen et al. H2003|l 
who report outward motions on both sides of the gap with 
similar velocities around 0.6 to 0.7 mas yr~^ corresponding 
to - 0.25 c. 

Besides the ionized (free-free absorbing) gas compo- 
nent there are multiple pieces of evidence for atomic and 
molecular gas in the central region of NGC 1052. H2O 
maser emission occurs towards the base of the western 
jet (Claussen et al. 1998) within the same region that 
is heavily affected by free-free absorption. Atomic hydro- 
gen is known to exist in NGC 1052 on various scales. Van 
Gorkom et al. |1986| imaged the distribution of the Hi 
gas with a resolution of 1" using the VLA. They report 
a structure three times the size of the optical galaxy. 
Recent VLBI observations resolved Hi absorption features 
towards the nuclear jet (Vermeulen et al. I2()()3|l . Finally, 
an OH absorption line was also detected by Omar et al. 
H2002|l and Vermeulen et al. I|2003|l but the distribution 
of the OH gas on parsec-scales has not been investigated 
so far. 

In this work we analyse the multi-frequency structure 
of NGC 1052 on parsec-scales between 5 GHz and 43 GHz 
in both total and linearly polarized intensity. In Sect. |21 
we describe briefly the observations with the Very Long 
Baseline Array (VLBA)^ and the data reduction. Section^ 
and Sect. E describe the modeling of the source structure 
with Gaussian components and the process of aligning 
the images at the four frequencies. The VLBA images of 
NGC 1052 themselves in total and polarized intensity are 
presented in Sect.[S| We analyse the frequency dependence 
of the observed core position in both jets in Sect. Eland 
Sect.ddiscusses the spectral analysis. In Sect.|Hlwe present 
the brightness temperature distribution along both jets of 
NGC 1052 and summarize our results in Sect.|^ 

2. Observations and data reduction 

NGC 1052 was observed on December 28th, 1998 with 
the VLBA at four frequencies (5 GHz, 8.4 GHz, 22 GHz, 
and 43 GHz) in dual polarization mode. The data were 
recorded with a bit rate of 128 Mbps at 2-bit sampling 

^ Napier 1199411 : the VLBA is operated by the National Radio 
Astronomy Observatory (NRAO), a facility of the National 
Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 



providing a bandwidth of 16 MHz per polarization hand 
(divided in two blocks of 16 0.5 MHz channels each). The 
total integration time on NGC 1052 was about one hour at 
5 GHz, 8.4 GHz, and 22 GHz each, and about six hours at 
43 GHz, to compensate the lower array sensitivity and the 
lower source flux density. 3C345 and 4C 28.07 were used 
as calibrators during the observation. The correlation of 
the data was done at the Array Operations Center of the 
VLBA in Socorro, NM, USA, with an averaging time of 
two seconds. All antennas of the array yielded good data, 
except the Owens Valley antenna, which did not record 
data at 8.4 GHz and 22 GHz. The data at 22 GHz and 
43 GHz suffered from some snow in the Brewster dish and 
rain at St. Croix. 

The data calibration and imaging were performed ap- 
plying standard methods using the programs AXVS and 
DIFMAP (Shepherd et al. 1997). The a-priori data cali- 
bration and fringe fitting were performed in AIVS us- 
ing the nominal gain curves measured for each antenna. 
Instrumental phase offsets and gradients were corrected 
using the phase-cal signals injected into the data stream 
during the data recording process. The data were aver- 
aged over frequency and exported from AIVS. Then, the 
data were read into difmap, edited, phase- and amplitude 
self-calibrated, and imaged by making use of the CLEAN 
algorithm. From a careful comparison of the uncalibrated 
and self-calibrated data, the absolute flux calibration can 
be (conservatively) estimated to be accurate on a level of 
< 10 percent. 

The data were corrected for the instrumental polar- 
ization of the VLBA using the method described by 
Leppanen et al. I|1995(l . The absolute values of the elec- 
tric vector position angles (EVPAs) at all four frequencies 
were calibrated using the source 3C 345, for which a large 
data base exists in the literature from which the EVPA 
in the core and jet regions of this source can be obtained 
(e.g., Ros et aL lM)n| . 

3. Modeling the source structure 

In VLBI imaging the absolute positional information is 
lost in the phase-calibration process. In the case of si- 
multaneous multi-frequency observations this means that 
a-priori it is not clear how the images at the different 
frequencies have to be aligned. The ideal method for over- 
coming this is to carry out phase-referencing observations, 
using a compact nearby object and using it to calculate the 
position of the target source relative to it (e.g., Ros 1200 j| . 
Our VLBA observations were not phase-referenced, so an- 
other way was used to register the four maps. We model 
fitted the visibility data with two-dimensional Gaussian 
functions. Those functions, called components, were cho- 
sen to be circular to reduce the number of free model 
parameters and to facilitate the comparison of the models 
at the different frequencies. The model fitting was per- 
formed in DIFMAP using the least-squares method. The 
errors were determined with the program erfit, a pro- 
gram from the Caltech VLBI data analysis package that 
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calculates the statistical confidence intervals of the fitted 
model parameters by varying each parameter. 

The fits were initially performed independently at the 
different frequencies to avoid biasing one by another. Once 
good fits for all four frequencies were obtained, a cross 
comparison of the resulting maps was made and the fits 
were modified to get a set of model fits as consistent as 
possible. Criteria for consistency were: 

— regions that show emission at adjacent frequencies 
should be represented by the same number of com- 
ponents, 

— extended components in the outer parts of the jets 
should become weaker at higher frequencies as they 
most likely represent optically thin synchrotron emit- 
ting regions, 

— the inner edges of both jets are expected to shift in- 
wards towards higher frequencies due to opacity ef- 
fects, i.e., the corresponding model components might 
have no low frequency counterpart, 

— optically-thin features should not show positional 
changes with frequency. 

Table^gives the parameters of the final models for the 
four frequencies. The most distant component in the east- 
ern jet was labeled as A 1, the adjacent inner one as A 2, 
and so on. The western jet was divided into three parts 
following the convention introduced above. The model fit 
component in the innermost part of the western jet were 
labeled as B2b, B2a, B 1 from east to west. Further out, 
the components C3b, C3a, C2, CI, and D follow. 

4. Image alignment 

The image alignment was performed in two steps, i) The 
fitted components at the different frequencies were cross- 
identified and the relative shifts between the models were 
determined, assuming frequency-independent positions of 
optically thin features, ii) An origin was determined from 
which absolute distances could be measured. This method 
of alignment is not a-priori definite. It is based on the as- 
sumption that optically-thin components have frequency- 
independent positions and that the cross identification of 
the model components is correct. However, all possible al- 
ternative identifications could be ruled out for consistency 
reasons (by making use of the component flux densities at 
adjacent frequencies). 

We used the two most distant components (A 1 and 
D) to align the 5 GHz and the 8.4 GHz models by as- 
suming that the mid-point between both components was 
spatially coincident. Because these two outer components 
are optically thin and not detectable above 8.4 GHz, the 
position of component A 7, which is relatively strong at 
all three frequencies, was used to align the three high- 
frequency models relative to the origin determined from 
the alignment of the two low- frequency models. As a nat- 
ural choice, the most probable position of the true center 
of jet activity, namely the center between the components 



Table 1. Model fit parameters 





Flux 








Id 


density 


Radius" 


P.A. 


FWHM 




[mjy] 


[mas] 


n 


[mas] 


5 GHz. Y„ri = : 


112 






Al 


80±2 


14.56±0.03 


66.5±0.1 


2.14±0.03 


A2 


54±2 


12.01±0.04 


63.6±0.1 


1.87±0.08 


A3 


60±2 


9.61±0.05 


65.3±0.3 


1.63±0.04 


A4 


150±7 


7.46±0.05 


71.7±0.2 


0.99±0.02 


A5 


338±37 


6.0±0.2 


72.8±0.1 


0.93±0.25 


A6 


447±48 


5.1±0.1 


72.8±0.2 


0.83±0.19 


A7/8 450±28 


4.0±0.1 


71.0±0.2 


0.69±0.08 


A9 


397±30 


3.1±0.2 


71.1±0.3 


0.20±0.07 


AlO 


140±37 


2.3±0.2 


70.9±0.4 


0.46±0.05 


C3 


137±3 


3.84±0.03 


-114.2±0.1 


1.03±0.06 


C2 


59±3 


4.75±0.06 


-115.74±0.3 


0.83±0.09 


CI 


24±1 


7.19±0.03 


-119.6±0.4 


1.85±0.10 


D 


67±1 


11.81±0.01 


-114.7±0.1 


2.24±0.05 


H.J. GHz. v,..^ = 


= 1 7fi 






Al 


38±1 


14.70±0.03 


66.5±0.1 


2.00±0.03 


A2 


37±1 


12.04±0.04 


64.0±0.1 


2.34±0.08 


A3 


36±1 


9.30±0.05 


65.9±0.3 


1.62±0.04 


A4 


83±4 


7.34±0.05 


71.9±0.2 


0.98±0.02 


A5 


72±8 


6.4±0.2 


72.7±0.1 


0.59±0.25 


A6 


360±26 


5.5±0.1 


72.8±0.2 


0.81±0.12 


A7 


396±22 


4.4±0.1 


72.3±0.2 


0.74±0.06 


A8 


146±8 


3.8±0.1 


69.7±0.2 


0.28±0.06 


A9 


360±21 


3.2±0.2 


71.3±0.3 


0.34±0.05 


AlO 


316±47 


2.7±0.2 


69.2±0.4 


0.31±0.03 


All 


137±25 


1.9±0.2 


74.3±0.4 


0.34±0.07 


B2a/b 


56±4 


0.60±0.03 


-119.8±0.1 


0.36±0.02 


Bl 


29±2 


1.26±0.03 


-119.9±0.1 


0.31±0.02 


C3 


20±1 


3.75±0.03 


-114.2±0.1 


1.06±0.08 


C2 


77±4 


4.86±0.06 


-114.9±0.3 


0.90±0.09 


CI 


21±1 


7.03±0.03 


-119.1±0.4 


1.68±0.10 


D 


48±1 


11.95±0.01 


-114.7±0.1 


2.21±0.05 


22 GHz. v„,, = 


71 






A6 


138±5 


5.48±0.03 


72.4±0.2 


0.98±0.04 


A7 


140±7 


4.39±0.01 


72.2±0.2 


0.65±0.03 


AS 


99±5 


3.72±0.01 


68.1±0.2 


0.47±0.03 


A9 


139±6 


3.09±0.01 


72.2±0.2 


0.32±0.02 


AlO 


144±5 


2.70±0.01 


68.5±0.2 


0.27±0.01 


All 


44±8 


1.98±0.04 


73.0±0.8 


0.29±0.07 


A12 


86±15 


1.53±0.04 


72.1±0.4 


0.26±0.06 


A13 


54±14 


1.20±0.04 


71.1±0.9 


0.22±0.06 


B2b 


341±14 


0.52±0.02 


-123±1 


0.21±0.01 


B2a 


151±14 


0.72±0.05 


-115±2 


0.29±0.02 


Bl 


31±2 


1.38±0.02 


-132±1 


0.3±1.7 


C3a 


50±4 


3.34±0.02 


-116.4±0.4 


0.58±0.03 


C3b 


86±6 


4.18±0.07 


-116±1 


1.95±0.16 


i,9 GHz. V,.., = 


83 






A6 


40±2 


5.36±0.03 


72.8±0.2 


0.92±0.05 


A7 


44±2 


4.39±0.01 


72.2±0.2 


0.51±0.03 


A8 


39±2 


3.76±0.01 


68.7±0.1 


0.46±0.03 


A9 


37±1 


3.177±0.003 


71.8±0.1 


0.22±0.01 


AlO 


70±1 


2.772±0.004 


68.4±0.1 


0.33±0.01 


All 


19±1 


1.95±0.01 


71.6±0.5 


0.32±0.03 


A12 


57±2 


1.499±0.004 


70.2±0.2 


0.22±0.01 


A13 


43±2 


1.221±0.006 


70.8±0.2 


0.12±0.01 


A14 


54±1 


1.017±0.003 


71.9±0.2 


0.08±0.05 


A15 


19±1 


0.47±0.01 


61±2 


0.21±0.03 


B2b 


225±2 


0.467±0.002 


-119.2±0.3 


0.224±0.002 


B2a 


31±2 


0.779±0.007 


-120.9±0.6 


0.22±0.01 


C3b 


27±2 


3.33±0.05 


-116.5±0.8 


0.95±0.08 



The radius is measured from the center between A15 and 
B2b (see Section 4.6). 
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NGC1052 
28/12/1998 



5GHz 



8.4GHz 



22GHz » 



43GHz 




20 



10 -10 

Relative R.A. (milliarcsec) 



-20 



Fig. 1. The aligned VLB A images 
of NGC 1052 at 5 GHz, 8.4 GHz, 
22 GHz, and 43 GHz. Contours and 
beam sizes are given in Table El 



A 15 and B 2b, was used as the origin (compare Sect. 
The component positions in Table ^ are given relative to 
this reference point. 

5. The brightness distribution in total and 
polarized intensity 

5.1. Total intensity imaging 

Figurenshows the aligned uniformly weighted total inten- 
sity images of NGC 1052 at the four frequencies and Fig.|21 
shows the images that result from applying Gaussian 
model-fit components to the measured visibilities. Tabled] 
gives the image parameters of Fig. ^ 

The basic source structure is formed of two oppositely 
directed jets divided by an emission gap. The source can 
be divided into four regions: region A makes up the whole 
(more or less continuous) eastward-directed jet emission; 
region B is not visible at 5 GHz, but becomes the brightest 
and most compact feature at high frequencies; regions C 
and D are bright only at low frequencies and become faint 
and diffuse at high frequencies. The main source charac- 
teristics are summarized in Tabic |21 A sub-division can 
be performed based on the Gaussian model fits, found to 
represent the source structure (see Fig. Region A is 
sub-divided into the components A 1 to A 15, region B 
into the components Bl, B2a, and B 2b, and region C 
into the components CI, C 2, and C 3. 

Whilst the eastern jet is only slightly curved at dis- 
tances from the gap larger than 4 mas the counterjet ex- 
hibits strong curvature most pronounced at 22 GHz. The 



Table 2. Total intensity image parameters 



V 


beam 




'S'peak 




rms 


Lowest Contour 


[GHz] 


[mas X mas 


■ °] 


[Jy/beam] 


[Jy] 


[m,Iy/bea 


Ti] [mjy/beam] 


5 


3.30x1.31, - 


3.74 


0.660*^ 


2.41 


0.26 


1.5 


8.4 


1.98x0.81, - 


3.87 


0.538*^ 


2.39 


0.25 


1.5 


22 


0.86x0.32, - 


7.63 


0.339"^ 


1.51 


1.20 


2.5 


43 


0.45x0.16, - 


7.93 


0.126^^ 


0.67 


0.67 


2.5 



Total flux density recovered in the image 
Corresponds to the A component 
Corresponds to the B component 



jets appear nearly symmetric in the 5 GHz and 8.4 GHz 
images, becoming asymmetric at higher frequencies. While 
the emission gap is most prominent at 5 GHz the images 
at 8.4 GHz, 22 GHz, and 43 GHz reveal jet components oc- 
cupying this gap region, leaving a smaller but still promi- 
nent emission gap. The frequency dependence of the VLBI 
"core" positions, the points were the eastern and the west- 
ern jet become optically thin can clearly be seen in Fig. ^ 
This so-called "core shift" will be analysed in detail in 
Sect. El 

5.2. Linearly polarized intensity imaging 

Figure |21 shows the linearly polarized intensity and its 
EVPA overlaid on the total intensity image at 5 GHz. A 
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5 




D 5GHz" 



8.4GHz 



22GHz 



43GHz 



A15 



■10 



Right Ascension (mas) 



Fig. 2. The aligned model-fit maps. 
The circles represent the full width 
at half maximum (FWHM) of the 
circular Gaussians. The vertical 
and oblique lines join the associ- 
ated components between different 
frequencies. The components A 1 
and D between 5 GHz and 8.4 GHz 
and A 7 between 8.4 GHz, 22 GHz 
and 43 GHz (joined by a thicker 
trace) were used for the align- 
ment. The lowest contours are 
1.9mJy/beam for the 5 GHz model, 
2.6mJy/beam for the 8.4 GHz 
model, 5.0mJy/beam for the 
22 GHz model and 3.5mJy/beam 
for the 43 GHz model. 



Table 3. Source characteristics on mas-scales 



Eastern Western 





Mean 


Jet 


Jet 


Gap 




V 


P.A. 


Length 


Length 


Width 


Notes 


[GHz] 




[mas] 


[mas] 


[mas] 




5 


65° 


30 


30 


~1.3 


A dominates; B totally 












absorbed 


8.4 


65° 


14 


14 


<0.5 


A dominates; B occupies 












the 5 GHz gap region 


22 


67° 


13 


6 


~0.5 


B dominates 


43 


70° 


5 


3.5 


~0.2 


B dominates; western jet 



weak 



region of linearly polarized emission is visible at the base 
of the eastern jet. The peak is about 3 mJy per beam and 
the EVPA is about 70°, roughly parallel to the jet. To 
decide whether the region of linearly polarized emission 
is resolved, slices along the jet axis of both the total in- 
tensity image and the polarization image were produced, 
using the task slice in ALVS (inlayed panel in Fig. 0) . In 
this plot, the polarized emission peaks about Imas offset 
from the total-intensity maximum. A comparison to Fig.|21 



ascribes this part of the jet to component A 10, which is 
optically thick at 5 GHz. The asymmetry of the polarized 
emission slice suggests the presence of at least two com- 
ponents. The polarized emission is thus slightly resolved 
and originates in the optically thick part of the eastern jet 
at 5 GHz. 

No linearly polarized emission at a flux density level 
above 1 mJy could be detected at the other three frequen- 
cies Usually, the degree of polarization rises with frequency 
since beam depolarization reduces the degree of polariza- 
tion at lower frequencies more strongly than at higher fre- 
quencies. Thus, one expects the eastern jet to exhibit po- 
larized emission from the region around component A 10 
and All at 8.4 GHz also. Assuming that the polarization 
has the same (flat) spectrum as the total intensity in this 
region, it should reach ^ 3 mJy per beam at 8.4 GHz (com- 
parable to the polarized flux at 5 GHz). However, this is 
not observed. At higher frequencies these components be- 
come optically thin (compare Sect. [7)) and thus fall in total 
intensity. Were their emission polarized to the same degree 
as at 5 GHz it would be < 1 mJy per beam and would 
thus lie below the detection threshold. The non-detection 
of linear polarization at the high frequencies in our obser- 
vations is consistent with the results of Middelberg et al. 
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Fig. 3. 5 GHz image of the polarized 
emission in NGC 1052 (EVPA su- 
perimposed on the total intensity). 
The inlaycd panel shows the total in- 



MiiMARC SEC tensity and polarized emission pro- 

files along the slice marked with the 
dashed line. 



H2()()4|l who report unpolarizcd emission from NGC 1052 
at 15 GHz down to a limit of 0.4%. 

6. Identifying the center of activity 

The symmetry between the jet and the counterjet con- 
strains the position of the central engine in NGC 1052. 
The cores of both jets are located at the distances to 
the central engine, Tc, where the optical depth r has 
fallen to '--^ 1. In a conical jet geometry this distance is 
given by: Vc ^ (see, e.g., Lobanov 11998(1 where 

kr = {{2a — 3)b — 2n — 2)/(5 — 2a), with a being the spec- 
tral index and b and n the power indices of the magnetic 
field and the density of the emitting particles: B ^ r^, 
N ^ r"^ fLobanov 11998(1 . Taking logarithms leads to: 

log(rc) = --^ log(z^) -I- const. (1) 

Measuring Tc at two frequencies allows one to determine 
kr in the corresponding region of the jet. For a freely ex- 
panding jet in equipartition (Blandford & Konigl il979i | ) . 
kr=\. The value of kr is larger in regions with steep pres- 
sure gradients and may reach 2.5, for moderate values of m 
and n (Lobanov .1998;) . If external absorption determines 
the apparent core position, comparable density gradients 
of the external medium can alter kr to values above 2.5. 

The values of kr deduced depend crucially on the abso- 
lute values of on the two sides and therefore on the as- 
sumed position of the central engine. Four scenarios have 
been tested with different reference points (see Fig. 0J). 
Table 0] gives the derived values of kr for each scenario. 
In each case the position of component A 14 has been as- 
sumed as the core of the eastern jet, rather than A 15. 
The latter is comparably weak and, thus, most likely does 
not represent the true jet core. A 15 might rather repre- 
sent a bright but heavily self-absorbed new jet component. 



Formally, however, the kr values for the scenarios 1,2, 
and 4 between 22 GHz and 43 GHz change to 1.4 ± 0.06, 
0.6 ± 0.03, and 1.0 ± 0.04 if A 15 instead of A 14 is used. 
For scenario 3, A 15 cannot be associated to the eastern 
jet core (compare Fig.0J. 

The area between the model components A 15 and B 2b 
is the most likely location of the central engine and the 
center between both components is a natural choice for 
its exact position (scenario 1). Shifting the reference point 
eastwards (scenarios 2 and 3) alters the values of kr into 
unphysically large regimes (requiring density gradients oc 
r"^*^ and higher). Assuming the true center of activity to 
be located more westwards (closer to B 2b, scenario 4), the 
values of kr derived are still acceptable. 

The positions of the bases of both jets at the differ- 
ent frequencies for the first case (scenario 1) are shown in 
Fig. The eastern jet has rather high values of kr be- 
low 22 GHz, although still in agreement with steep pres- 
sure gradients in the jet environment. Above 22 GHz kr is 
3.9±0.8, which is a good indicator for free-free absorption 
affecting the jet opacity. The western jet has values of kr 
as high as 6.8±2.7 between 22 GHz and 43 GHz, suggest- 
ing a large contribution from free-free absorption. 

The results from the core shift analysis support the 
picture of a free-free absorbing torus covering mainly the 
inner part of the western jet and also a smaller fraction of 
the eastern jet. The true center of activity in NGC 1052 
can be determined to lie between the model components 
A 15 and B 2b, with an uncertainty of only ~ 0.03 pc. 

7. Spectral analysis 

Spectral information can be derived from multi-frequency 
VLBI data in two ways. The first approach is to use the 
knowledge of the proper alignment of the total intensity 
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Fig. 4. Image with the model fitting results of the 43 GHz 
data. The innermost jet components are labeled. The pu- 
tative locations of the central engine are indicated for four 
different scenarios (see discussion in the text). 



Eastern Jet 




3 5 10 20 50 

Frequency [GHz] 

Fig. 5. Core positions in the two jets at the different fre- 
quencies for scenario 1 (see Fig. 0J . Table 0] provides the 
values of kr for the other two scenarios. The inset panel 
shows the core locations in the source at the four frequen- 
cies. 



images to derive images of the spectral index between two 
adjacent frequencies. The second approach is to derive 
spectra of the model fit components. The approaches are 
somewhat complementary. The latter gives a handy num- 
ber of component spectra, which can be analysed in detail, 
whereas the spectral index imaging gives the full course 
of the spectral index along the jet axis. Particularly, this 
yields spectral information at parts of the jet that are not 
adequately represented by Gaussian model components. 
The results of both approaches will be presented in this 
section. 



Table 4. Values of kr for the four different putative cen- 
ters of activity 



(5-8.4) GHz 


(8.4-22) GHz 


(22-43) GHz 


Scenario kr, esist ^r,west 


^r,east ^r,west 




1 2.7±1.9 - 


2.1±0.5 6.8±3.0 


4.1±0.8 6.2±2.2 


2 2.4±1.7 - 


1.8±0.4 9.8±4.3 


3.1±0.6 9.3±3.4 


3 2.1±1.5 - 


1.4±0.3 12.9±5.2 


2.2±0.4 12.5±4.6 


4 3.0±2.1 - 


2.5±0.6 3.7±1.6 


4.9±1.0 3.0±1.0 


Table 5. Spectral index image parameters 


1/1 U2 


Common Beam 


Cutoff 


[GHz] [GHz] 


[mas X mas, °] 


[mjy/beam] 


5 8.4 


2.64x1.06, -3.81 


2 


8.4 22 


1.42x0.57, -5.75 


1 


22 43 


0.66x0.24, -7.78 


1 



7.1. Spectral index imaging 

Spectral index images have been produced using the 
AXPS task COMB^. For this, the total intensity has been 
reimaged with appropriate tapering of the (u, w)-data to 
match the resolutions at adjacent frequencies. Information 
below ImJy/beam was discarded in both input images. 
Table El gives the restoring beams and the other relevant 
parameters of the derived images. The spectral index im- 
ages are shown in Fig. El 

The main feature in each image of the spectral index is 
an optically thick inner edge of the innermost part of both 
jets. Outwards along the jets, the spectral index tends to 
decrease. The spectral index a exceeds the value of 2.5, 
the theoretical upper limit for synchrotron self-absorption, 
on both sides of the gap in the (5-8.4) GHz spectral in- 
dex image and on the eastern side of the gap also between 
8.4 GHz and 22 GHz. High values of a are also reached 
on the western side of the gap in the (8.4-22) GHz spec- 
tral index image and on both sides between 22 GHz and 
43 GHz, which however do not exceed the critical value of 
2.5. Towards the outer jet regions the spectral index typ- 
ically falls to values < — 1 in the whole frequency range 
between 5 and 43 GHz. This suggests that the emission 
from these outer jet regions far from the gap is optically 
thin above 5 GHz. 

7.2. Spectral analysis of the model fit components 

The flux densities of the model fit components (see 
Figure O can be used to derive spectra of i) the whole 
parsec-scale structure, ii) the two jets separately and 

At each pair of pixels of the two (coinciding) images, the 
spectral index a {S oc is calculated from the flux densities 
per beam Si and 5*2 at the two frequencies i^i and 1^2: a = 
(log Si - log S2)/(log Ui - log U2) 
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Fig. 6. Spectral index images of NGC 1052. The contours 
and beam parameters are given in Table |5| Note the dif- 
ferent scales in the three images and that the images are 
centered on component B2. 



Ill) the model fit components themselves. In Fig. [7| the 
total spectrum of the parsec-scale structure of NGC 1052 
between 5 GHz and 43 GHz is shown, as well as the spectra 
of both jets separately. All values were obtained by adding 
up flux densities^ of model components (A 1 to A 15 for 
the eastern jet and B 2b to D for the western one). 

At the high frequencies both jets show similar spectra 
with an optically thin decrease of the flux density above 
22 GHz. The spectral index in this regime is around —1 in 
both jets, with the eastern jet being significantly stronger 
than the western one in agreement with the interpretation 
that the eastern jet approaches the observer whereas the 
western jet is the counterjet. However, below 22 GHz the 
spectra of the two jets differs substantially. The eastern jet 
spectrum remains optically thin above 8.4 GHz and flat- 
tens around 5 GHz. The western jet, on the other hand, 
exhibits a sharp decrease in flux below 22 GHz. Although 
the spectral index limit of 2.5 is not exceeded, synchrotron 



* We estimated the statistical standard uncertainty of all 
quantities derived from the model component flux densities. 
The calculated errors are shown in the plots, when they are 
not negligible. 




• 


Total Flux 


* 


Eastern Jet 




Western Jet 
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Frequency [GHz] 

Fig. 7. Spectrum of the parsec-scale structure of 
NGC 1052 between 5 GHz and 43 GHz (total and for both 
jets separately). Errorbars are shown only if they exceed 
the symbol size. 

self- absorption seems very unlikely to be responsible for 
the turnover of the spectrum because of two reasons. First, 
the similarity of the spectra of eastern and western jet at 
high frequencies suggests similar intrinsic physical prop- 
erties on both sides so that the self absorption frequency 
should not differ by a factor of three. Second, the kine- 
matical analysis of Vermeulen et al. H2003|l shows the jet 
axis to be close to the plane of the sky and the motions 
to be only weakly relativistic. There is little evidence for 
strong Doppler boosting. Moreover, such an effect should 
shift the turnover frequency of the counterjet to lower fre- 
quencies, rather than to higher ones. In Sect. we will 
use the moderate velocity differences, on the one hand, 
and the more pronounced brightness-temperature differ- 
ence of jet and counterjet, on the other hand, to constrain 
the orientation of the jet-counterjet system. 

If the spectrum of the western jet is decomposed into 
the spectra of the different jet regions B, C and D (see 
Fig. n and Fig.|SJ) it turns out that the region B, the in- 
nermost region on the western side, is responsible for the 
inversion of the western jet spectrum below 22 GHz. Its 
spectrum between 5 GHz and 8.4 GHz is highly inverted 
with a spectral index > 4 since it is not detectable at all at 
5 GHz. Such a spectral index cannot be due to synchrotron 
self absorption but indicates a region of external absorp- 
tion towards the core of the western, receding jet. The 
highly inverted spectrum of this component in NGC 1052 
and the necessity of an external absorber was first men- 
tioned by Kellermann et al. (,L999.) and later confirmed by 
Kameno et al. H2001fl and Vermeulen et al. H2003fl . 

It is not clear whether the eastern jet is also affected by 
external absorption. Kameno et al. (|2001,1 proposed a ge- 
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Fig. 8. Spectra of the four jet regions A, B, C and D (see 
Fig. n for the definition of the jet regions). The arrows 
indicate that component C is partly resolved at 43 GHz 
and component B is not detected at 5 GHz. Thus, the 
shown values represent only a lower and upper limit , 
respectively. Errorbars are shown only if they exceed the 
symbol size. 

ometry in which the obscuring torus, covers 0.7 mas of the 
western jet and 0.1 mas of the eastern jet. Emission from 
A 11 is not detected at 5 GHz. Estimating (conservatively) 
a flux density of lOmJy at 5 GHz yields a spectral index 
of ~ 5. This would suggest that the innermost part of the 
eastern jet is also strongly affected by external absorption. 
However, if the flux densities of the two components A 10 
and All are summed at 5 GHz and 8.4 GHz, one obtains 
a spectral index of 2.3, which is close to that expected 
for pure synchrotron self-absorption. Consequently, due 
to the small separations between the inner components 
in the eastern jet, we cannot finally judge from this ap- 
proach whether free-free absorption plays an important 
role. Kameno et al. 2003 determined a spectral index of 
3.28 ±0.27 for the base of the eastern jet between 1.6 GHz 
and 4.8 GHz and attribute this to free-free absorption in 
the obscuring torus. 

The jet-to-countcrjet ratio of NGC 1052 can be deter- 
mined from the model fitted flux densities of both jets. 
Figure El shows the ratio of the flux densities of the model 
components on either side of the gap as a function of fre- 
quency. At high frequencies the jet-to-counterjet ratio is 
~ 1.5 and starts rising towards lower frequencies below 
22 GHz. At 8.4 GHz the jet is brighter than the counter- 
jet by a factor of ^ 5 and at 5 GHz the jet outshines the 
counterjet even by a factor of ~ 7. The jet-to-counterjet 
ratio is known to keep rising at lower frequencies up to 
50 at 2.3 GHz (Kameno et al. IMIUl . This suggests that 
a much bigger part of the western jet is covered by the 
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Fig. 9. The jet-to-counterjet ratio determined from the 
flux densities of the model fit components as a function 
of frequency. Errorbars are shown only if they exceed the 
symbol size. 

absorber than apparent from our high-frequency data. 
Alternatively, curvature effects may play an important 
role. In this case, the counterjet bends away from the ob- 
server at a distance of a few milliarcseconds from the core 
and its radiation is Doppler de-boosted. 

8. Brightness-temperature gradients 

In this section the brightness-temperature gradients along 
both jets in NGC 1052 are discussed. Starting from the 
assumption that the magnetic field B, the electron density 
iV, and the jet diameter D can be described by power-laws: 

B (xr'';N (xr";D (xr"^ (2) 

it can be derived that the brightness temperature, Tb, is 
expected to fall with increasing distance from the jet base 
like: 

Tbtxr^ , (3) 

with s < 0. For this, one assumes optically thin syn- 
chrotron emission with an emissivity ji, oc ncB{v / vb)'^ 
(Krolik ' 19991 Sect. 9.2.1) and, implicitly, a constant 
Lorentz factor of the emitting electrons. Under these as- 
sumptions the brightness temperature distribution is de- 
termined by the jet-geometry (d = 1 for a conical jet, 
d < for a coUimated jet, and d > for a decelerat- 
ing jet), the course of the magnetic field and the particle 
density via 

s = d + n + h{l-a) , (4) 

where a is the spectral index. Typically, straight and con- 
tinuous jets exhibit values of s ~ —2.5 on parsec-scales 
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(Kadler et al., in prep). The search for deviations from 
this power-law dependence provides a tool to find regions 
in a jet which are affected by external absorption or by 
abrupt changes of the jet parameters. 

For a non-thermal source, the brightness temperature 
is a frequency-dependent quantity (see e.g., Condon et al. 

with Si, the flux density of the source, v the observing 
frequency, and Q the apparent diameter at half maxi- 
mum. In principle, the observing frequency has to be cor- 
rected for the source redshift, but due to the small dis- 
tance of NGC 1052 (z = 0.0049) this correction is negli- 
gible. According to eq. brightness temperatures have 
been computed for the model fit components at all four 
frequencies. The uncertainties in Tb were computed using 
Gaussian error propagation from the errors in flux density 
S and FWHM of the model components 8. We assumed 
conservative values of 10% for AS* and 20% for A6. 




1 10 
Distance Along Eastern Jet [mas] 



Fig. 10. Brightness temperature distribution along the 
eastern jet of NGC 1052. Above the cut-off distance of 
~ 2.5 mas the data can be fitted by a power-law (see 
Table |H1 for the fit parameters) with a power-law index 
of ~ —4. The offset between the four data sets of bright- 
ness temperatures reflects the frequency dependence of Tb 
(see equation EJ. 



8.1. The Tb distribution along the eastern jet 

Figure [TUl shows the brightness temperatures in the east- 
ern jet as a function of distance from the central engine. 
The latter was assumed to lie at the center between the 
components A 15 and B 2b (see Sect. El for a detailed dis- 
cussion of the central engine position in NGC 1052). Tb 
rises towards the center following roughly a r~^-law (com- 
pare Table EJ . Between 2 mas and 3 mas from the central 
engine, however, there is an abrupt decrease of Tb (dashed 
line in Fig. I10|l . This "cut-off" is present at all four fre- 
quencies although the inner components (only visible at 
22 GHz and 43 GHz) exhibit again a rise in Tb. The value 
of the innermost component A 15 falls significantly below 
the extrapolation of the curve defined by A 12, A 13, and 
A 14. 



Table 6. Tb gradient along the eastern jet 



v 

[GHz] 


[K] 


s 


5 

8.4 

22 
43 


2.1+li X w 

5.3tii X 10'' 
1.9lJ-^ X 10" 


-4.1 ±0.4 
-4.1 ±0.3 
-3.8 ±0.3 
-4.1 ±0.8 



Note: We performed a power-law fit to Tb beyond a distance 
from the central engine r = 2.5 mas as Tb(r) — Ti, x r" , where 
is the extrapolated brightness temperature at 1 mas and s 
is the power index. 



8.1.1. What causes the cut-off of the Tb-distribution? 

We discuss three possible origins of the frequency- 
independent cut-off of Tb in the eastern jet: 



Free-free absorption: At a distance of ~ 2.5 mas from the 
central engine the overhanging edge of the obscuring torus 
might start to obscure a substantial fraction of the jet and 
thus reduce the brightness temperature of its components 
via free-free absorption. If this interpretation was correct 
the pronounced frequency dependence of the optical depth 
due to free-free absorption (ry oc i/~^'^) should be mea- 
surable: the free-free absorbed (i.e., observed) flux density 
Su.abs depends on the intrinsic flux density S,^, the optical 
depth of the absorber tj and the observing frequency v as 



5*1^, abs — Si/ 



(6) 



The calculated values of Tf are given in TableElfrom which 
it is obvious that no pronounced effect of decreasing opac- 
ity with increasing frequency is present. This makes the 
interpretation as free-free absorption unlikely. 

Synchrotron self-absorption: The synchrotron emission 
from a freely expanding, relativistic jet is self-absorbed 
at distances smaller than 



jet 



c + 1 



1 



obs 



l/{n+e6+l) 



(7) 



due to the change of the optical depth along the jet 
(Lobanov[ 



(8) 



n+eb+l re 



M. Kadler et al.: The Twin-Jet System in NGC 1052 



11 



6 



1 






▲ 5 GHz 

□ 8.4 GHz 
• 22 GHz 

43 GHz 


















1 


T 


: 


B2b B2a 1 Bl 


^ 

^3 C2 


Cl 




1 






10 


Distance 


Along Western 


Jet [mas] 



Fig. 11. Brightness temperature distribution along the 
western jet of NGC 1052. No simple power-law can ap- 
proximate the data, probably due to strong curvature ef- 
fects. 

Here, e = 3/2 — a, i/ is the observing frequency, 0obs 
is the observed opening angle of the jet (thus, implicitly 
assuming d = 1; compare Equation EJ, if jet is a bootstrap 
constant describing the jet conditions at a certain char- 
acteristic distance (Lobanov 1998), and 6j is the jet bulk 
Doppler factor. Acceleration and deceleration of the flow 
may affect the dependence of Equation Q , and ultimately 
even cause the observed abrupt decrease of the brightness 
temperature at shorter distances from the central engine 

(Fig. Da. 

The effect of deceleration can be accounted for by 
assuming (5j oc r-'^ , with / < 0, and ^obs oc r'^, with c > 0. 
With these assumptions, the length of the self-absorbed 
portion of the jet becomes oc j^^, with 



^ e{b + f) + n + c+l' ^ ' 

The peak-positions of the brightness temperature dis- 
tribution in the eastern jet fFig. I10|l imply ^ > —1/10 ( 
note that C < 0) for the shallowest possible slope of r^,. 
For a typical set of assumption about the jet flow (b — —2, 
n = —1, a = —1) this requires / < —8 — 2/5c, which is 

Table 7. Optical depth values - eastern jet 



Id 5 GHz 8.4 GHz 22 GHz 43 GHz 



AlO 3.2 ±1.1 

\ — 3.0 ±0.9 2.5 ±1.0 3.0 ±1.9 

A12 — — 2.6 ±0.9 2.3 ±2.3 

A13 — — 3.6 ±1.1 2.2 ±2.7 

A14 _ _ _ 1.9 ±3.3 

A15 _ _ _ 8.2 ±6.9 



Note: Values for a distance from the central engine smaller 
than r = 2.5 mas 



implausible as it implies extremely strong deceleration 
{Sj oc r^^ and higher) of the flow. This scenario can there- 
fore be ruled out. 

Steep pressure gradients: Frequency-independent local 
maxima of the brightness-temperature distribution along 
the jet can also be produced by strong density and pres- 
sure gradients at the outer edge of the nuclear torus and, 
to a lesser degree, by the gradients of the magnetic field 
strength (Lobanov 1998). Such gradients may result in 
rapid changes of synchrotron self-absorption and exter- 
nal free-free absorption of the jet emission, both of them 
increasing the opacity at shorter distances from the nu- 
cleus. These two factors together can in principle explain 
values of C > ^1/10, if density gradients at the outer 
edge of the torus are stronger than n « —10 {N oc r"). 
For the expected size of the nuclear torus rtor ~ 0.2 pc to 
0.3 pc (corresponding to the distance of 2 mas to 3 mas at 
which the brightness- temperature decrease sets on), and 
typical densities in the absorbing torus (tih ^ lO^cm"'^ 
to lO^^cm"'^, Cassidy & Raine 1993) and in the nuclear 
medium of the host galaxy (tt-h ~ 10^ cm''^ to lO^cm"^, 
Ferguson et al. 11997(1 . the density increase should occur 
on scales of ~ 2.5rtor- This seems to be the most plausi- 
ble explanation for the observed brightness temperature 
changes in the eastern jet. 

Standing shocks: The approach to describe the jet param- 
eters B, N, and D with simple power laws (Equation IJ)) 
assumes, implicitly, a quasi-stationary, continuous jet flow 
without perturbations. This model cannot describe the in- 
fluence of shocks in the jet flow on the brightness tempera- 
ture distribution. Qualitatively, the occurence of a cut-off 
in the Tb-distribution at the location of a region of en- 
hanced linearly polarised emission (compare Sect l5.2|l is 
in agreement with the expected observational signature of 
a standing shock. Part of the bulk flow energy of the jet 
can be converted into magnetic energy of the jet plasma 
in such a standing-shock region. This might cause Kjct to 
increase abruptly. Beyond the shock, the jet flow might be 
quasi-stationary and thus well described by our model. 

8.1.2. Derivation of the spectral index 

The frequency dependence of the brightness temperatures 
in the eastern jet of NGC 1052 (see Table enables 
us to derive directly the value of the spectral index a: 
Equation [S] shows that the brightness temperature at a 
given distance from the central engine depends on the 
spectral index as Tb oc i/"^"*"", if optically thin synchrotron 
emission is assumed. A linear regression of the brightness 
temperature at a distance of Imas {T{i'); see TableEl as 
a function of frequency v yielded T(:/) oc i/-3 0±o.i_ xhus, 
the spectral index is a = — l.OiO.l. This is in good agree- 
ment with the results of the spectral index imaging (see 
Sect. EH). 
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8.1.3. The gradients of particle density and magnetic 
field 

The fitted sizes of the Gaussian components at the four 
frequencies do not show significant deviations from a con- 
ical jet-structure so that d = 1 can be assumed. Applying 
this to Equation leads to 

n + 2b 5 . (10) 

In the most simple scenario of a conical expanding jet 
with a well defined particle energy distribution function 
(i.e. without cooling; n = —2) and equipartition between 
magnetic energy and particle energy (6 = n/2), a relation 
n + 2h = —4 is expected. Our result thus implies that at 
least one of these two assumptions does not hold strictly. 
If one assumes that equipartion holds, energy losses (e.g. 
due to adiabatic expansion) might steepen the particle en- 
ergy distribution effectively to n —2.5. An alternative 
explanation would be that the magnetic flux of a domi- 
nant longitudinal component of the magnetic field density 
along the jet axis {ttD^Bz) might be conserved so that 
b = —2. In the extreme case of no transverse magnetic 
field component this would lead to n + 2b —6. 

In principle, the detection of polarized emission in the 
eastern jet (see Sect. 15 .211 allows the direction of the dom- 
inant component of the magnetic field to be determined. 
The linearly polarized emission originates in the region 
around component A 10, thus at the edge of the most 
compact component of the obscuring torus. Faraday ro- 
tation in this region is expected to be large so that we 
cannot directly derive a dominant transverse magnetic 
field from the alignment of the EVPAs with the jet axis. 
Polarimetric observations at multiple closely separated 
frequencies around 5 GHz would be necessary to derive the 
amount of Faraday rotation, which in combination with 
the measured opacity due to free-free absorption Tf ad- 
ditionally would allow to disentangle the particle density 
and the magnetic field (the Faraday rotation at a given 
frequency is proportional to the particle density rie, while 
Tf oc nl). 

8.2. The Tb distribution along the western Jet 

The brightness- temperature distribution along the west- 
ern jet of NGC 1052 is shown in Fig. ^2 Here, the situation 
is more complex than on the eastern side. The western jet 
is less continuous than the eastern jet, which might be due 
to stronger curvature, particularly, between the regions C 
and D (see the counterjet structure at 22 GHz in Fig. QJ. 
Since the equations derived above assume a straight jet 
geometry without abrupt bends, we do not try to approx- 
imate the brightness-temperature distribution by a single 
power-law. However, the course of the brightness tempera- 
ture along the western jet tells something about the orien- 
tation of the jets and the structure of the obscuring torus 
towards the receding jet. 

Tb decreases rapidly outwards between 3 and 12 mas. 
At a distance of about 7 mas, however, in the area of 



component CI, the brightness temperature has a local 
minimum at 5 GHz and 8.4 GHz, which might be due to 
Doppler de-boosting due to a bend away from the line of 
sight. Inwards of 3 mas, in region B, where the strongest 
effects of free-free absorption were detected (see Sect, d 
the Tb distribution is flatter and the peak brightness tem- 
peratures do not exceed values of a few times 10^° K. 

Component B2 is strongly absorbed and, thus, 
its brightness temperature is substantially reduced at 
8.4 GHz, dropping even below the peak value at 22 GHz 
(as expected due to its inverted spectrum discussed 
above). This can be used to estimate the absorbing column 
density towards component B2, the area of the strongest 
effects of free-free absorption in NGC 1052. Assuming an 
intrinsic symmetry between jet and counterjet we can use 
the measured frequency-dependence of T'(z^) from the east- 
ern jet of NGC 1052 to estimate the intrinsic brightness ra- 
tio of the counterjet at 8.4 GHz and 22 GHz. Considering 
the slightly different fitted slopes (compare Table ^ we 
derive from the best sampled region between 3 mas and 
8 mas along the eastern jet that the brightness tempera- 
ture at 8.4 GHz should exceed the value of Tb at 22 GHz by 
a factor of 15 to 20. Component B 2 has approximately the 
same brightness temperature at 8.4 GHz and at 22 GHz 
(mean value of B 2a and B 2b) which means that its flux 
density at 8.4 GHz is reduced by a factor of ~ e~^'^-e~^ °, 
i.e., t| *^~^^(B2) - 2.7 to 3.0. This is in agreement with 
the results of Kameno et al. ('2001'). These authors derive 
an optical depth of 300 at 1 GHz in the region correspond- 
ing to our component B2, yielding Tf ^ 3.4 at 8.4 GHz. 
The optical depth due to free-free absorption is given by 
(e.g., Lobanov .l998jl 

= 30- 10i6xr-i-35z/-2iTi2 . (11) 

For r/ ~ 3 at 1/ = 8.4 GHz at a temperature of T = 10^ K, 
and a length of the absorber of i = 0.3 pc (comparable to 
the extent of the absorbing region in the plane of the sky) 
we derive a density of n = 2.4 • 10^ cm~^ and an absorbing 
column density of 2.2 • 10^^ cm~^ towards component B 2. 
Depending on the unknown ionization fraction of the torus 
materiaP this value is consistent with various X-ray ob- 
servations of NGC 1052, which imply a (model-dependent) 
column density of iVn ~ lO^^cm"^ to A^h = lO^'^cm"^ 
towards the unresolved nuclear X-ray core (Weaver et al. 
nM?l Guainazzi et al.HnMI Kadler et al.|2I3|). Because 
of the smaller absolute difference between the correspond- 
ing values of Tb at 5 GHz and 8.4 GHz and the relatively 
larger uncertainties we do not derive opacities for the outer 
components of the western jet. However, from the inspec- 
tion of Fig. ^2 it is clear that the influence of free-free 
absorption is weak beyond 4 mas west of the nucleus. 



^ Free-free absorption is an indicator of the column den- 
sity A'^o of free electrons in an ionised medium, rather than 
the column density of neutral hydrogen A^h- In a fully ionised 
medium, A''e = Ah holds. 
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8.3. Orientation of the jet-counter jet system 5. 

The ratio of brightness temperatures in the jet and in the 
counterjet at the same distances from the central engine 
can be used to constrain the angle to the line of sight of 
the jet/counterjet axis: 

= ("i^V"" , (12) 

Sj_ ^ 1 + /3 cos 6> 

6cj l-/?cos6i ' ^ 6. 

where Sj and 6cj are the Doppler factor of the jet and the 
counterjet, respectively. The mean ratio of the measured 
values of Tb on the western side and the fitted value at 
the corresponding distance on the eastern side (calculated 
for components between 3 mas and 5 mas distance from 
the center where free-free absorption effects are expected y 
to be small) at all four frequencies is 9 ± 2. Assuming a 
spectral index of —1, this gives (3cos9 = 0.35lg;Q4. Since 
/3 = 1 is an upper limit for the jet speed this results in a 
maximum allowed angle of ~ 72°. The minimum allowed 
angle derived by Vermeulen et al. ( 2003 ) is ^ 57° for which 
/3 = 0.64. Thus, the jets are constrained to lie at an angle g 
to the Hne of sight between 57° and 72°. 

9. Discussion and Summary 

1. VLBI imaging of NGC 1052 exhibits a parsec-scale 
"twin-jet" structure, matching the standard model of 
AGNs if the two jets are oriented close to the plane of 9. 
the sky. 

2. We present accurately aligned, high-quality VLBI im- 
ages of NGC1052 at 5 GHz, 8.4 GHz, 22 GHz, and 
43 GHz, the associated spectral index images between 
the adjacent frequencies and spectra of the various jet 
regions and model fit components. 

3. The core of the western jet has a highly inverted 
spectrum with a spectral index well above 2.5, the 
theoretical upper limit for synchrotron self absorp- 
tion, which was first mentioned by Kellermann et al. 
I|1999|) and later confirmed by Kameno et al. ( 200T) 
and Vermeulen et al. ( 2003). Qualitatively, we confirm 
the results of those authors, particularly, the increas- 
ing opacity in the inner tens of milliarcseconds of the 
eastern jet. 

4. We analysed the frequency dependence of the observed 
VLBI core position in both jets and found another 
clear signature of free-free absorption at the core of the 10. 
western jet. The shift rate with frequency is too low 

to be explained in terms of synchrotron self absorp- 
tion alone, while the core shift rate on the eastern side 
can still be explained under the assumption of steep 
pressure gradients increasing the synchrotron opacity. 
From the determination of these core shift rates we 
obtain an independent measurement of the position 11. 
of the central engine of NGC 1052 with an accuracy as 
high as ~ 0.03 pc, superior to a kinematical derivation. 



We find a sharp cut-off of the brightness temperature 
distribution along the eastern jet. Neither synchrotron 
self absorption nor free- free absorption can explain this 
behaviour. The most plausible explanation is that we 
see an effect of steep pressure gradients at a transition 
regime between the external pressure-dominated jet 
regime and a more or less freely expanding jet regime. 
Thus, the sharp cut-off of the brightness temperature 
distribution marks an observational sign of the over- 
hanging edge of the obscuring torus. 
We used the ratio of the observed brightness temper- 
atures in the jet and the counterjet to constrain the 
angle to the line of sight of the "twin-jet" system. 
Together with the information from the kinematical 
study of Vermeulen et al. (j2[)()3|l . the angle to the line 
of sight can be determined to lie between ~ 57° and 
- 83°. 

The spectral index of the synchrotron jet emission 
was found to be —1. This result comes independently 
from the frequency dependence of the brightness- 
temperature distribution in the eastern jet and from 
the imaged spectral index at large distances from the 
core. 

Either equipartition between the magnetic energy 
and the particle energy or the assumption of a sin- 
gle, well defined particle energy distribution with- 
out cooling is violated in the parsec-scale eastern jet 
of NGC 1052. Alternatively, a conserved longitudinal 
component might dominate the magnetic field on these 
scales. 

We find a region of linearly polarized emission at the 
base of the eastern jet. The EVPA of the polarized 
emission cannot be evaluated directly since Faraday 
rotation at the center of this galaxy is expected to be 
large. We find no linear polarization of the source at 
higher frequencies. The simplest explanation for this 
behavior is that the different layers of the jet have a 
different degree of polarization and that the emission 
at higher frequencies originates in an inner, unpolar- 
ized layer of the jet. The higher column density at 
the more centrally located jet base at higher frequen- 
cies could cause higher Faraday depolarization than it 
does at 5 GHz. This idea is supported by the fact that 
the jet base at 5 GHz coincides with the edge of the 
obscuring torus. This scenario can be tested with ob- 
servations at longer wavelengths and especially around 
5 GHz, where the depolarization is expected to start 
dominating. 

Most likely, a combination of free-free absorption, syn- 
chrotron self absorption, and the presence of steep 
pressure gradients determine the parsec-scale radio 
properties of NGC 1052. An analytical model fit to the 
observed spectrum of only one single absorption mech- 
anism model seems not to be satisfactory to represent 
the true physical situation. 

The absorbing column density derived from the de- 
graded brightness temperature of the western jet core 
is ~ 2 X 10^^ cm~^, in good agreement with the value 
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obtained from X-ray spectroscopy (Kadler et al.[ 
This suggests that the nuclear X-ray emission (which 
is unresolved for all X-ray observatories currently in 
orbit) originates on the same scales which are imaged 
by VLBI and underlines the importance of combined 
future radio and X-ray observations of NGC 1052. 

Acknowledgements. We are grateful to T. Beckert, A. Kraus, 
T. P. Krichbaum, and A. Roy for many helpful discussions and 
suggestions. We thank the referee, Seiji Kameno, for his careful 
reading of the manuscript and his suggestions, which improved 
the paper. M. K. was supported for this research through a 
stipend from the International Max Planck Research School 
(IMPRS) for Radio and Infrared Astronomy at the University 
of Bonn. 



References 

Blandford, R. D. & Konigl, A. 1979, ApJ, 232, 34 

Cassidy, I. & Raine, D. J. 1993, MNRAS, 260, 385 

Claussen, M. J., Diamond, P. J., Braatz, J. A., Wilson, A. S., & 

Henkel, C. 1998, ApJ, 500, L129 
Condon, J. J., Condon, M. A., Gisler, C, & Puschell, J. J. 

1982, ApJ, 252, 102 
Ferguson, J. W., Korista, K. T., & Ferland, G. J. 1997, ApJS, 

110, 287 

Guainazzi, M., & AntoneUi, L. A. 1999, MNRAS, 304, L15 
Kadler, M., Kerp, J., Ros, E., et al. 2004, A&A, 420, 467 
Kameno, S., Sawada-Satoh, S., Inoue, M., Shen, Z., & Wajima, 

K. 2001, PASJ, 53, 169 
Kameno, S., Inoue, M., Wajima, K., Sawada-Satoh, S., & Shen, 

Z. 2003, Publ. of the Astron. Soc. of Australia, 20, 134 
Kellermann, K. I., Vermeulen, R. C, Cohen, M. H., & Zensus, 

J. A. 1999, BAAS, 31, 856 
Knapp, G. R., Faber, S. M., & Gallagher, J. S. 1978, Astron. 

J., 83, 139 

Krolik, J. H. 1999, Active galactic nuclei : from the central 
black hole to the galactic environment, Princeton, NJ: 
Princeton University Press 

Leppanen, K. J., Zensus, J. A., & Diamond, P. J. 1995, AJ, 
110, 2479 

Lobanov, A. P. 1998, A&A, 330, 79 

Middelberg, E., Roy, A. L., Bach, U., et al. 2004, in 
Future Directions in High Resolution Astronomy, ASP 
Conf. Ser., J.D. Romney & M.J. Reid (eds.), in press, 
(astro-ph/0309385) 

Napier, P. J. 1994, lAU Symp. 158: Very High Angular 
Resolution Imaging, J. G. Robertson & W. J. Tango (eds.), 
p.117 

Omar, A., Anantharamaiah, K. R., Rupen, M., & Rigby, J. 

2002, A&A, 381, L29 
Ros, E., Zensus, J. A., & Lobanov, A. P. 2000, A&A, 354, 55 
Ros, E. 2004, in Future Directions in High Resolution 

Astronomy, ASP Conf. Ser., J.D. Romney & M.J. Reid 

(eds.), in press, ( astro-ph/0308265 1 
Shepherd, M. C. 1997, ASP Conf. Ser. 125: Astronomical Data 

Analysis Software and Systems VI, 6, 77 
van Gorkom, J. H., Knapp, G. R., Raimond, E., Faber, S. M., & 

Gallagher, J. S. 1986, AJ, 91, 791 
Vermeulen, R. C, Ros, E., Kellermann, K. I., et al. 2003, A&A, 

401, 113 

Weaver, K. A., Wilson, A. S., Henkel, C, & Braatz, J. A. 1999, 

ApJ, 520, 130 
Wrobel, J. M. 1984, ApJ, 284, 531 



